On the possibility of an excitonic insulator at the semiconductor-semimetal transition 
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We calculate the critical temperature below which an excitonic insulator exists at the pressure- 
induced semiconductor-semimetal transition. Our approach is based on an effective-mass model 
for valence and conduction band electrons interacting via a statically screened Coulomb potential. 
Assuming pressure to control the energy gap, we derive, in the spirit of a BEC-BCS crossover 
scenario, a set of equations which determines, as a function of the energy gap (pressure), the 
chemical potentials for the two bands, the screening wave number, and the critical temperature. 
We (i) show that in leading order the chemical potentials are not affected by the exciton states, (ii) 
verify that on the strong coupling (semiconductor) side the critical temperatures obtained from the 
linearized gap equation coincide with the transition temperatures for BEC of non-interacting bosons, 
(iii) demonstrate that mass asymmetry strongly suppresses BCS-type pairing, and (iv) discuss in 
the context of our theory recent experimental claims for exciton condensation in TmSeo.45Teo.55. 

PACS numbers: 71.35.-y, 71.30.+h, 71.28.+d, 75.20.Hr, 75.30.Mb 



I. INTRODUCTION 

More then four decades ago Mott realized that a 
semimetal (SM) with a low density of free charge carriers 
is unstable against an insulating, semiconducting (SC) 
state, because the Coulomb interaction binds conduction 
band electrons and valence band holes to excitons. At 
around the same time, Knox Q made the important ob- 
servation that a SC whose energy gap is smaller then 
the exciton binding energy has to be unstable against 
the spontanous formation of excitons. Understood as a 
global instability, the pairing of electrons and holes is ex- 
pected to lead to a macroscopic, phase coherent quantum 
state or condensate - the excitonic insulator (EI) - sepa- 
rating below a critical temperature the SC from the SM. 
Theoretically 0, 0, S B S H ES E3 , the instability 
should occur in any material which can be tuned from an 
indirect gap semiconductor to a semimetal, for instance, 
by applying pressure, uniform stress, or by optical pump- 
ing. Yet, experimental efforts to establish the EI in real 
compounds largely failed. 

It is only until recently that detailed studies of the 
pressure- induced SC-SM transition in TmSeo.45Teo.55 by 
Wachter and coworkers strongly suggested the existence 
of an EI [H [H [H El 

In particular, the anomalous 
increase of the electrical resistivity in a narrow pressure 
range around 8 kbar indicated the appearance of a new 
phase below 250 K ^| . Wachter and coworkers sug- 
gested that this new phase might be an "excitonic insu- 
lator", and, assuming pressure to modify only the energy 
gap Eg, constructed a phase diagram for TmSeo.45Teo.55 
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in the E g — T plane [13J . Later they found in the same 
material a linear increase of the thermal diffusivity below 
20 K and related this to a superfluid exciton state 

E3 

Both excitonic phases are located on the SC side of the 
SC-SM transition in TmSeo.45Teo.55. 

In the present paper we perform the first theoretical 



analysis of these astonishing experimental findings, stay- 
ing strictly within the bounds of the concept of an EI 
as a superfluid condensate of electron-hole pairs, simi- 
lar to a condensate of Cooper pairs in a superconduc- 
tor. (We exclusively use the acronym EI to denote this 
state and refer to condensation when electron-hole pairs 
enter this state.) In contrast to the early theoretical 
works H H H i, H H H EH E3 , we (i) self-consistently 
calculate the phase boundary of the EI, (ii) analyze the EI 
in the spirit of a BEC-BCS crossover, and (iii) investigate 
also the halo of the EI, that is, the region surrounding 
the EI on the SC side. 

The creation of an exciton condensate is usually at- 
tempted by optical pumping of suitable semiconductor 
structures So far, without success. The main ob- 

stacle is the far-off-equilibrium situation caused by op- 
tical excitation. In contrast, pressure-induced genera- 
tion of excitons occurs at thermal equilibrium, which 
is much more favorable for condensation. Semiconduct- 
ing, pressure-sensitive mixed valence materials, such as 
TmSeo.45Teo.55, offer therefore a very promising route to- 
wards exciton condensation. It is even conceivable to use 
this class of materials for implementing recent propos- 
als of coherent transport across EI-SM junctions |17lll8| . 
The validation of Wachter and coworkers' experiments 
would thus have a tremendous impact on the field of ex- 
citon condensation. 

The excitonic instability is driven by the Coulomb in- 
teraction leading to pairing of conduction band (CB) 
electrons with valence band (VB) holes. Of particu- 
lar importance is therefore to self-consistently determine 
its weakening when the external parameter, for instance 
pressure, pushes the material from the SC side, with 
only a few thermally excited charge carriers available for 
screening, to the SM side, with a huge number of charge 
carriers. The strength of the Coulomb interaction deter- 
mines on which side of the SC-SM transition the system 
is. As long as the Coulomb interaction is strong enough 
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FIG. 1: Qualitative sketch of the phase boundary of an EI 
with equal band masses. (Adopted from Ref. 



to support excitons, that is, as long as the exciton bind- 
ing energy is positive, the material is on the SC side. The 
vanishing of the binding energy (Mott effect Q) defines 
then the SC-SM transition. 

Depending from which side of the SC-SM transition 
the EI is approached (see Fig. the EI typifies either a 
Bose-Einstein (BE) condensate of tightly bound excitons 
(SC side) or a BCS condensate of loosely bound electron- 
hole pairs (SM side) 0,|2(J- A characteristic feature of 
BEC of excitons on the SC side is that formation of exci- 
tons and phase coherence (condensation) occur at differ- 
ent temperatures ji^, El ■ While excitons already exist 
below the temperature Tm set by the Mott effect, the EI, 
understood as a genuine condensate of excitons, occurs 
only below the critical temperature T c < Tm- Thus, on 
the SC side, the EI is embedded in a region containing ex- 
citons in addition to VB holes and CB electrons. We call 
this region halo. On the SM side, in contrast, Cooper- 
type electron-hole pairs form and condense at the same 
temperature T c . The difference between an EI sitting, 
respectively, on the SC and SM side of the SC-SM tran- 
sition has important consequences for the interpretation 
of the TmSeo.45Teo.55 data. 

After the seminal studies by Leggett Comte and 
Nozieres [i^, and Nozieres and Schmitt-Rink j^H, the 
transition from BE to BCS condensation has been ex- 
tensively studied in the past, mostly with an eye on 
short-coherence length superconductors [i^, HE 0, H^, 
Efil E^R I2II but also electron-hole gases in semiconduc- 
tors 12(1 13 (1 . and, more recently, trapped atomic Fermi 
gases |3lL l32j | have been analysed from this point of view. 
From these analytical and numerical investigations it is 
known that the transition from BE to BCS condensa- 
tion is smooth. It is also known that diagrammatic ap- 
proaches usually capture the crossover only when the 
BCS equation for the order parameter is augmented by 
an equation for the chemical potential. At zero tem- 
perature, it suffices to calculate both quantities in mean- 
field approximation [l9ll20|. At finite temperatures, how- 
ever, the chemical potential has to be at least determined 
within a T-matrix approximation which accounts for the 
excitation of collective modes, i.e. in the context of an 
EI, of excitons with finite center-of-mass momentum |2l| . 
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In order to discuss the EI at the pressure-induced SC- 
SM transition in terms of a BEC-BCS crossover, we need 
therefore not only equations for the order parameter and 
the screening wave number (to account for the Mott ef- 
fect) but also for the CB electron and VB hole chem- 
ical potentials. Because pressure controls the chemical 
potentials, which in turn are constrained by charge neu- 
trality, the meanfield approximation turns out, even for 
finite temperatures, to be sufficient on the (weak cou- 
pling) SM side and the (strong coupling) SC side, in con- 
trast to what one would expect from related diagram- 
matic studies of the crossover problem in high-T^ mate- 
rials and atomic Fermi gases E3, 0, EM 
Using the Thouless criterion |3jj (see also 
mine the transition temperature for BEC directly from 
the normal phase electron-hole T-matrix, we verify that 
the transition temperatures obtained from the meanfield 
equations indeed coincide on the SC side with the tran- 
sition temperatures for BEC of non-interacting bosonic 
excitons. From the normal phase electron-hole T-matrix 
we furthermore extract the temperature Tm above which 
excitons cease to exist and construct a mass-action-law 
to determine the composition of the El's halo. 

As far as Wachter and coworkers' experiments [lil H3I 
IliL HrI | are concerned we come to the following conclu- 
sions. The phase boundary constructed from the resistiv- 
ity data does not embrace the EI. Instead, it most proba- 
bly reflects that part of the halo of the EI where excitons 
prevail over free electrons and holes and give rise to an 
efficient additional scattering channel which leads to the 
observed resistivity anomaly. More remarkable is how- 
ever that from our theoretical results we would expect 
the EI - that is, the macroscopic, phase coherent conden- 
sate - exactly in the temperature range where the linear 
increase of the thermal diffusivity was observed 0|. In 
analogy to liquid He 4, Wachter and coworkers suggest 
that the increase could be due to the second sound of a 
superfluid exciton liquid. We, on the other hand, find the 
EI entirely on the SC side, where, within our approxima- 
tions, it constitutes a BEC of ideal bosonic excitons. The 
BCS side of the EI, where it is unclear whether superflu- 
idity can be realized or not 0, is strongly suppressed 
for finite temperatures because of the large asymmetry 
between the CB and VB masses in TmSeo.45Teo.55. Of 
course, an ideal gas of bosons cannot become superfluid 
and thus cannot feature a second sound. In reality, how- 
ever, excitons interact and could, when the interaction is 
dominantly repulsive, give rise to it. Thus, Wachter and 
coworkers [l5| may have seen an exciton condensate. 

In the next section we first introduce an effective-mass 
model on which our calculation is based and then give 
a description of our calculational scheme. In section 
III we present numerical results for the phase bound- 
ary T c (E g ) of the EI, first for equal band masses, where 
we obtain a steeple-like T c (E g ), directly reflecting the 
different condensation mechanisms on the SC and SM 
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side, and then for asymmetric band masses, applicable 
to Tm[Se,Te] compounds, where we find a strongly sup- 
pressed T c (E g ) on the SM side and thus an EI sitting 
almost entirely on the SC side of the SC-SM transition. 
We then relate in section IV our results to the experi- 
mental data for the pressure induced SC-SM transition 
in TmSeo.45Teo.55 and conclude in section V. 

II. FORMALISM 
A. Model 

The excitonic instability arises because of the Coulomb 
attraction between electrons in the lowest CB (i = 1) 
and holes in the highest VB (i = 2), with an indirect 
energy gap separating the two bands. Since the spin 
algebra would unnecessarily mask the many-body theo- 
retical concepts we want to discuss, we focus on spinless 
fcrmions. Only in section III, where we make contact 
with experiments, we include the spin in our calculation. 

Keeping only the dominant term of the Coulomb in- 
teraction, which we assume to be statically screened, an 
effective-mass model for studying the (spinless) Wannier- 
type EI is 

H = J2 ei(k)ct kCi;k + \ Wp(<lV(-q). (!) 

k,i q 

with p(q) = £V k cj k+q Ci,k the total charge density and 
V s (q) = (Aire 2 / eo) / (q 2 + q 2 ). The screening wave num- 
ber q s depends on the CB electron and VB hole density 
and will be determined self-consistently in the course of 
the calculation; eo is the background dielectric constant. 
The momenta k for CB electrons and VB holes are mea- 
sured from the respective extrema of the bands which are 
separated by half of a reciprocal lattice vector w = K/2. 
We consider the case where both the VB and the CB 
band have one extremum. Assuming isotropic effective 
masses rrii, the band dispersions read (Ti = 1 throughout 
the paper) 

ei(k) = E g + £l {k)-n, e 2 (k) = -e 2 (k)- M -£ (k), (2) 

with £i(k) = k 2 /2m.;, E g the bare energy gap which can 
be varied continuously through zero under pressure, and 
fi the chemical potential. The band structure (0) refers to 
the unexcited crystal at T = with an empty CB and a 
full VB whose selfenergy Eo (k) has to be subtracted [13 • 
For model to be applicable for the description of 
an EI we have to assume that through the SC-SM tran- 
sition all parameters, except the energy gap E g and the 
screening wave number q Sl vary weakly and can be kept 
constant. Since model treats an indirect gap semi- 
conductor as a direct gap semiconductor, effects due to 
the indirect gap |icj . in particular, multi- valley effects 



(when the valleys are not identical) and the balancing of 
the finite momentum w of the electron-hole pairs by the 
lattice leading, for instance, to (small) density modula- 
tions in the EI phase are beyond the scope of the present 
paper. 

B. Meanfield approximation 

In formal analogy to the theory of superconductivity, 
we employ matrix propagators (in the band indices i = 
1,2) 

Gy(k,T) = -(T T c k>i (r)c k> .(0)), (3) 

with < r < f3 and T T [...] the time-ordering operator 
with respect to r. The diagonal elements Gn denote the 
(normal) propagators for electrons in the CB and VB, re- 
spectively, whereas the off-diagonal (anomalous) propa- 
gators describe the phase coherence of electron- hole pairs. 
Finite anomalous propagators signal the condensate and 
thus the EI. 

In Matsubara space, the matrix propagator satisfies a 
Dyson equation, 

G(k)=g(k)+g(k)Z(k)G(k), (4) 

with a diagonal matrix 

g i3 (k) = [iuj n - e i (k]~ 1 5i j , (5) 

denoting the bare propagator and a 2x2 selfenergy ma- 
trix S, which we assume to be given by a skeleton expan- 
sion, that is, in terms of the fully dressed propagator G. 
The variable k = (k, iu> n ) denotes a four- vector with k a 
momentum and uj n = (2n+I)7r//3 (n integer) a fermionic 
Matsubara frequency. 

Separating the selfenergy into diagonal (normal) and 
off-diagonal (anomalous) parts and suppressing the k de- 
pendence, we rewrite the Dyson equation as 

g = g + g^ A G, (6) 
g = 9 + gi: N g (7) 

with Yj A and Y> N the anomalous and normal selfenergy, 
respectively. 

The separation of the Dyson equation is particularly 
useful for the calculation of the transition temperature 
T c . We first note from Eq. J7J that g is diagonal. Thus, 
Eq. (jfJJl reduces to 

r -i = ( Qxx - fe|Sf 2 | 2 [SnG^r 1 \ 

' " I [feGn4] _1 g^-g^\ 2 )- 1 > 

In the vicinity of the transition temperature T c the 
anomalous selfenergy becomes vanishingly small and 
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(a) 



(b) 



+ 



(c) 



FIG. 2: (a) exchange selfenergy E x to be used in the mean- 
field approximation, (b) ladder selfenergy E L additionally 
used in the T-matrix approximation, and (c) ladder equa- 
tion for the T-matrix. Solid (dashed) lines denote matrix 
propagators (statically screened Coulomb potentials). 



Eq. (jHJ) can be linearized with respect to E A . The anoma- 
lous propagator G\i is then given by 



G12 — Q\\Q22^ 



22^12! 



(9) 



where the propagators Qn satisfy Eq. (JJJ) with E w cal- 
culated from Qu alone; G12 and G21 do not enter E^ 
anymore. Thus, within the linearized theory, the calcu- 
lation of the normal propagators is decoupled from the 
calculation of the anomalous propagators. From Eq. 
we see moreover that G12 7^ is equivalent to E^ 2 ^ 0. 
Accordingly, the anomalous selfenergy can be considered 
as an order parameter. 

To proceed we have to specify the selfenergies E" and 
Ti A . In the meanfield approximation both are given by 
the exchange energy H x shown in Fig. Applying 
standard diagrammatic rules |3fij . we find 



E£(k) = 



dk' 



(2tt) 3 

for the normal selfenergies and 



y s (k-k> F (e,(k') (10) 



A(k) 



dk' 

WY 3 



V s (k-k')Q(k')A(k') 



with 



Q(k) = 



n F (e 2 (k)) - n F (ei(k)) 
ei(k) - e 2 (k) 



(11) 



(12) 



for the anomalous selfenergy A(k) = E£j(k), where n F 
denotes the Fermi function. 

Note, the normal selfenergies contain no Coulomb-hole 
term because we replaced directly in the Hamiltonian the 
bare Coulomb interaction by the statically screened one 
(static screening approximation); the correlation energy 
is thus partly suppressed. As a consequence, excitonic 
phases are energetically more favorable than electron- 
hole liquid phases j33|, independent of the mass ratio 



and the number of valleys, and in contrast to, for in- 
stance, what one finds in the quasi-static approxima- 
tion |23, which performs the replacement (with a mod- 
ified screened Coulomb potential) after the selfenergies 
have been calculated in the random-phase approxima- 
tion [3^. Neither approximation, however, should be 
used to study the competition between the two phases, 
because both do not give a reliable correlation energy 
for intermediate densities, that is, close to the SC-SM 
transition, where excitons should be included in the cor- 
relation energy ^0]. In addition, scattering processes 
beyond our model, for instance, between valleys or with 
phonons, also affect the competition between excitonic 
and electron-hole liquid phases. An investigation of the 
competition between the two is thus a complex task, be- 
yond the scope of the present paper, which simply as- 
sumes excitons to be favored from the start. The static 
screening model is thus sufficient. 

The single particle energies entering Eqs. 1)1 0|l and l|12(l 
are the self-consistent solutions of 



e i (k) = e i (k)+E^(k) 



(13) 



The main effect of the exchange energy on the particle 
dispersions is a rigid, k-independent energy shift which 
we incorporate into chemical potentials for the CB elec- 
trons and VB holes. Writing the renormalized disper- 



e i (k) = (-l) i + 1 ( £i (k)- Mi ) (14) 
yields for the electron and hole chemical potentials 

Mi = V- E g- Ae i> (I 5 ) 
/i 2 = Ae 2 - (U, (16) 

with energy shifts satisfying 



Ae, ; 



dk 



V fl (k)n F (£i(k) -fa). (17) 



To derive Eq. H17[l for i = 2, the selfenergy Eo of the full 
VB has to be calculated within the same approximation 
and then drops out. Together with Eq. ltT7|) . Eqs. (TSJ 
and (|16|) effectively define an electron-hole representation 
which is convenient for the calculation of the screening 
parameter q s . 

In order to obtain a closed set of equations for the 
three unknown parameters /1, fix, and /i 2 , we augment 
Eqs. lfl"5l) and l(TB|) by the condition of charge neutrality 
which forces CB electron and VB hole densities, given in 
the meanfield approximation, respectively, by 



ni 



n 2 



dk 

(2^)3 

dk 

(2^)3 



n F (ei(k) - fix), 
n F (e 2 (k) - (Jq), 



(18) 
(19) 
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to be equal. Thus, charge neutrality leads to the con- 
straint 



[n F (ei(k) - m) - n F (e 2 (k) - /z 2 )] = 0, (20) 



which we use to eliminate [i in Eqs. (|15f) and l|16(l . For 
that purpose, we add Eqs. ((T5|) and l(TB|) 



Ml + M2 = --Eg + Y1 J (^3 y ^ k ^ F 



( £i (k)-/i i )(21) 



and combine this equation with the charge neutrality con- 
straint (|2(J|) . Physically, the sum of the CB electron and 
VB hole chemical potentials is the negative of the renor- 
malized energy gap: 



E g = ei(0)-e 2 (0) 



-Mi - M2- 



(22) 



With the individual chemical potentials for each 
species at our disposal, the screening parameter q s is 
given by 



47re 2 



d 



d 



i Q — m + q — n 2 

e \Ofli 0^2 



(23) 



Once Mi;M2, an d <?s are calculated, the kernel of the 
gap equation (fTTJ) is known in the whole E g — T plane. 
For a given energy gap E g , we can thus determine the 
transition temperature T c as the highest temperature for 
which the homogeneous integral equation (|llf) has a non- 
trivial solution A(k) ^ 0. 



C. T- matrix approximation 

The meanfield theory presented in the previous sub- 
section determines T c (E g ) on the strong coupling (SC) as 
well as on the weak coupling (SM) side. The transition 
temperatures on the SC side turn out to coincide with the 
BEC transition temperatures for a non-interac ting boson 
gas of excitons (see section III), as it should |t3.l2lj. From 
the perspective of diagrammatic BEC-BCS crossover the- 
ories UllUmiHliSlillllllii this is however sur- 
prising because the meanfield theory does not account for 
excitons existing on the SC side already above the tran- 
sition temperature T c (E g ). These excitons are expected 
to at least modify the chemical potentials However, 
the charge neutrality constraint, inherent in a pressure- 
driven BEC-BCS crossover, leads to a cancellation of the 
leading order corrections to the chemical potentials. 

Diagram matic BEC-BCS crossover theories 0,113,0, 
Ell27ll2i.l3lll32]| usually leave the gap equation (JT]} un- 
changed but calculate the charge densities from normal 
propagators dressed not only with the exchange selfen- 
ergy E x (Fig. |2K) but also with the ladder selfenergy E L 
(Fig. |2t>), which takes, via the electron- hole T- matrix 



(Fig. Et), normal phase excitons into account. Expand- 
ing the spectral functions to which the total selfene rgy 
gives rise to with respect to the imaginary part of E L |41| 
(see appendix A), the CB electron (VB hole) density can 
be written as n\ — n{ + n\ (n 2 — + n\) with the 
superscripts 'f and 'b' denoting, respectively, the part 
of the total density coming from free electrons (holes), 
with a modified dispersion but a spectral weight unity, 
and the part arising from bound electron-hole pairs. 
The charge neutrality condition H2(J[> becomes therefore 
n{ + n\ = t%2 + n\- Now it is important to realize (or 
to verify by calculation, see appendix A) that whenever 
a CB electron participates in an exciton, a VB hole has 

to do the same. Thus, n\ = n^, and the charge neutral- 

f — f 

ity condition reduces to n{ = n 2 which, except for the 
modifications of the dispersions due to E L has the same 
form as in the meanfield approximation. 

The above argument is independent from the approx- 
imation used to obtain the electron-hole T-matrix. It 
is only based on charge neutrality and the fact that 
an exciton has to contain both an electron and a hole. 
The calculation shows that it is crucial to employ the 
skeleton expansion. As a result, the single particle dis- 
persions are modified by E L and thus by excitons (see 
Eq. HA.45I) - ). The ladder approximation (Fig. Et) we use 
to calculate the T-matrix does however not consistently 
account for feedback effects of this kind; for that purpose, 
additional diagrams |26| or a systematic mode-coupling 
approach [2^, |2j| should be implemented. We consider 
therefore feedback of excitons to single particle states 
as higher order effects and ignore all of them, includ- 
ing the one already present in the ladder approximation. 
The charge neutrality condition in the ladder approxi- 
mation reduces then to the meanfield charge neutrality 
condition. Note, feedback effects induce interactions be- 
tween pairs jl^. Neglecting them implies therefore to 
stay within the framework of a non-interacting, ideal gas 
of excitons. How to extract from the skeleton expansion 
(and not from an algebraic mapping which is question- 
able in the vicinity of the SC-SM transition) the interac- 
tion between electron-hole pairs, which not only contains 
the repulsive core at short distances because of the Pauli 
exclusion but also the weakly attractive tail at large dis- 
tances, is to the best of our knowledge an unsolved prob- 
lem and beyond the scope of the present paper. 

In appendix A we calculate the charge densities within 
the ladder approximation shown in Fig. using the sep- 
arable approximation for the normal phase electron-hole 
T-matrix described in appendix B. Thereby, we proof by 
calculation that n\ — fii reduces to n{ = n 2 ■ Because the 
screening parameter q s , on the other hand, accounts only 
for screening due to free charge carriers it has to be calcu- 
lated from n{ and n 2 . Thus, even for finite temperatures, 
the meanfield approximation of the previous subsection 
is on par with T-matrix based BEC-BCS crossover theo- 
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ries |2Jl|2j|2j|25l|2Jl|2J,|3il|32j. 

As an additional check, we now deduce the transition 
temperature on the SC side directly from the normal 
phase electron-hole T-matrix (which takes excitons with 
finite center-of-mass momentum into account) using the 
Thouless criterion [33| (see also 0) and verify that the 



J 



T c (E g ) calculated within the meanfield approximation in- 
deed coincides with the transition temperatures obtained 
from the T-matrix. Thereby, we also verify that the ne- 
glect of the feedback of excitons on single particle states 
leads to a non-interacting bosonic gas of excitons. 



In the approximation depicted in Fig. [2J;, the normal phase electron-hole T-matrix is given by 

/dk" 
j^V s (k - k")Q(k"; q, ffi„)A 12 (k", k'; q, iQ r 



with an electron-hole pair propagator 

Q(k; q,ifl n ) 



1 - np(fi(k) - fxi) - n F (e 2 (k - q) - /jg) 
Eg +ei(k) +e 2 (k-q) -ifl n 



(24) 



(25) 



Here, we used Eq. (|14|l to express renormalized dispersions in terms of chemical potentials pi and definition 1221) 
to introduce the renormalized energy gap; Q„ — 2rnr/P (n integer) are bosonic Matsubara frequencies and q denotes 
the center-of-mass momentum of an electron- hole pair (and thus of an exciton). 

Taking advantage of the fact that the most important part of the normal phase electron-hole T-matrix originates 
form the exciton state, we calculate in appendix B the T-matrix within a separable approximation. As a result, we 
obtain 

2 

(26) 



A 12 (k, k'; q, = $(|k - ^±q|) • D x (q, itt n - E g - • g(\k< - ^q|) 



with a renormalized exciton propagator 

Dx(q_,iQ n ) 

and form factors 



-1 



?0 7l + S + M x (q,zfi„) 



7(P) 



-V s (p,p')x(p') = 



2m 



D 



x(p) 



(27) 



(28) 



defined in terms of a screened exciton wavefunction \ and a screened exciton binding energy B (see below and appendix 
B). Here, V s (p,p') denotes the screened Coulomb potential averaged over the angle between p and p', M = mi + m 2 
is the total mass, and m = mim 2 /M is the reduced mass of an electron-hole pair. 

The exciton propagator l|27|) describes an exciton in a medium consisting of a finite density of CB electrons and 
VB holes. The screened exciton binding energy B already accounts for the screening due to spectator particles not 
participating in the exciton. The spectators' Pauli blocking, on the other hand, is included through the selfenergy 



M x (q,ifi n ) 



dp lh+B\ 



■x\p) 



n F (ei(p+ ^-q) 



TO 2 

Mi) + n F (e 2 (p - — q) - p 2 ) 



(29) 



To avoid unnecessary numerical work, we replace in Eq. I|28|l the statically screened Coulomb potential by the 
Hulthen potential 



V s (r) 



e 2 exp (-q s r) 



Vh(t) 



2? s 



e exp (2q s r) - 1 ' 



(30) 



which is a good approximation as long as the screening wave number q s is not too large. Indeed, in the limit q s — > 0, 
the Hulthen potential reduces to the unscreened Coulomb potential. With the replacement (|30l) . Eq. (|28|) can be 
solved analytically [42^ |. The screened exciton wavefunction and the screened exciton binding energy for ax • q s < 1 
are, respectively, given by 



x(p) 

B 



ly/ira 3 x (l - (a x ■ q s ) 2 ) 



[(1 - a x ■ q s ) 2 + (a x • p) 2 ] ■ [(1 + a x ■ q s ) 2 + (a X ■ p) 2 ] 

Rx(l - ax ■ q s ) 2 , 



(31) 
(32) 
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where Rx = l/2ma| is the exciton Rydberg and 
ax = £o/me 2 is the exciton Bohr radius (e bare elec- 
tron charge). For ax ■ q s > 1, the Hulthen potential 
does not support a bound state and we have to set 
X = B = in that parameter range. Notice, in the 
limit q s — > 0, Eqs. (pTTJl and J32J), respectively, reduce to 
the wave function and binding energy of an isolated exci- 
ton (no medium effects) as obtained from Eq. (|28|l with 
V s replaced by the angle averaged unscreened Coulomb 
potential. 

The pole of the exciton propagator l|27|) determines 
the analytical structure of the T-matrix. Physically, it 
gives the exciton binding energy B(q) renormalized by 
screening and phase space filling. Assuming a weak q- 
dependence, B(q) ~ B(0) = B, with the renormalized 
exciton binding energy B determined from 

B = B + ReMx(0,ifl n -> -B + irf). (33) 

In terms of B the exciton propagator can be rewritten as 

D x (q.,iSl n ) = (34) 
i\ l n + B 

with an exciton spectral weight defined by 

Z x = l--^ReM x (0,n + ir))\ n= _ S . (35) 

The Thouless criterion states that BEC of excitons oc- 
curs when the normal phase electron-hole T-matrix di- 
verges at q = and id n — 0. Using Eq. I|26|l this implies 
D X X {Q, —Eg) = 0. With Eq. the transition temper- 
ature T c (E g ) is thus given by 

B{E g ,T c ) — Eg(Eg,T C ), (36) 

where we displayed the dependence of the renormalized 
binding energy and the renormalized energy gap on the 
bare energy gap and the temperature. 

In order to make the physical content of Eq. Ij36|l more 
transparent, we rewrite it in terms of chemical potentials. 
Recalling the definition (1221) of the renormalized energy 
gap, Eq. ffify becomes 

= Mi + = -B, (37) 

which is equivalent to the criterion for BEC of a non- 
interacting Bose gas: The transition occurs when the 
chemical potential y,x of the bosons reaches the bottom 
of the boson band leading to a macroscopic occupation 
of the q = state. 

The Thouless criterionper se jsjjj can be used to deter- 
mine T c {E g ) on the SC and the SM @ side. However, 
by construction, Eq. i|2"6*|) describes the normal phase 
electron-hole T-matrix only on the SC side, where ex- 
citons exist, and the separable approximation for the T- 
matrix is applicable. Thus, we can use the Thouless cri- 
terion only on the SC side. 



The exciton binding energy obtained from Eq. Ij33(l en- 
ables us also to determine the phase boundary Tm(E 9 ) 
between the SC and the SM and thus the halo of the 
EI, that is, the region between Tm(E 9 ) and T c (E g ) on 
the SC side, where excitons, CB electrons and VB holes 
coexist. As originally suggested by Mott .1], we use the 
vanishing of B (Mott effect) as a criterion for the SC-SM 
transition. Accordingly, Tm{E 9 ) is given by 

B(E g ,T M ) = 0. (38) 

III. RESULTS 

To construct the phase boundary of the EI, we dis- 
cretize the gap equation and determine, for a given 
energy gap E g , the temperature T c (E g ) for which the de- 
terminant of the coefficient matrix of the resulting linear 
set of equations vanishes. For each pair (E g ,T) we sup- 
ply the chemical potentials /ii and together with the 
screening wave number q s by finding the simultaneous 
roots of Eqs. (gDJ, igTJ, and To obtain the SC- 

SM phase boundary, we first determine, from Eq. 13311 . 
the renormalized exciton binding energy B as a function 
of Eg and T, again providing for each pair (E g ,T) the 
chemical potentials and the screening wave number. Af- 
ter that, we search for a given E g for the temperature 
Tn,i(Eg) which satisfies Eq. I|38|) . In a similar way, we 
determine on the SC side, as a cross-check of the results 
obtained from the linearized gap equation, the critical 
temperature T c (E g ) from the Thouless criterion (|36|l . 

We consider an isotropic system, angles can be thus 
integrated analytically. All integrals over the magnitude 
of the momentum are done by Gaussian integration ex- 
cept the |k|-integration appearing in the gap equation, 
which we perform by product integration, with the an- 
gle averaged Coulomb potential V s (k, k') as the weight 
function. The logarithmic singularity of the kernel at 
k = k' , which arises on the SC side for T — > 0, can 
thus be handled without problems. Recall, the separable 
approximation for the electron-hole T-matrix, on which 
the calculation of the renormalized exciton binding en- 
ergy B is based, uses the Hulthen potential instead of 
the statically screened Coulomb potential. All results are 
presented in exciton units, measuring energies and tem- 
peratures in exciton Rydbergs Rx and lengths in exciton 
Bohr radii ax- 

The phase boundary T c (E g ) for an EI with equal band 
masses (a = mi/ni2 = 1) is presented in Fig. The 
solid line shows T c (E g ) obtained from the linearized gap 
equation The solid circles on the SC side, on the 

other hand, give the T c (E g ) deduced from the Thouless 
criterion (|36|l for BEC. The two critical temperatures co- 
incide almost perfectly supporting our claim that even for 
finite temperatures the meanfield theory captures BEC 
on the SC side. The small deviations between the two 
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FIG. 3: Phase boundary of an EI with equal band masses 
(a = mi/mz = 1). The solid line (full circle) denotes the 
critical temperature obtained from the linearized gap equa- 
tion (Thouless criterion for BEC). The dotted line shows part 
of the phase boundary on the SC side when no screening due 
to thermally excited charge carriers is taken into account. The 
thick (thin) dashed line indicates the SC-SM transition, de- 
fined by the vanishing of the renormalized (screened) binding 
energy, taking screening and Pauli blocking (only screening) 
into account. 



curves as we come closer to the SC-SM transition (thick 
dashed line) originate from the discrepancy between the 
Hulthen and the screened Coulomb potential when the 
screening wave number increases. 

From Fig. |3 we see that above T * 0.12 the EI is un- 
stable and an ordinary SC-SM transition occurs (thick 
dashed line), here defined by the vanishing of the renor- 
malized exciton binding energy B. For T < 0.12, we find 
a steeple-like phase boundary, which strongly discrimi- 
nates between E g > and E g < 0. For E g > 0, T c (E g ) 
first increases rather rapidly with increasing E g , reaches 
a maximum, and then decreases to zero at E g = 1, 
the critical energy gap, above which the EI cannot ex- 
ist 0, H II @, LLJJ. For E g < 0, in contrast, T c {E g ) 
initially drops relatively fast, stays almost constant in a 
narrow _E 9 -range, and then decreases monotonously with 
further decreasing E g . Notice, in contrast to the qual- 
itative phase boundary of Fig. ^ the calculated T c (E g ) 
reaches the BCS asymptotics only at very large band 
overlaps, far away from the SC-SM transition. On the 
scale of the exciton binding energy, we find no simple 
exponential dependence of the transition temperature on 
the band overlap \E g \. 

The steeple-like shape of the phase boundary reflects 
the different character of the EI when it is approached 
from the SM and the SC side, respectively. Deep on the 
SM side (E g < and \E g \ > 1, not shown in Fig. OJ), 
the EI constitutes a BCS condensate of loosely bound 
electron-hole pairs whose small binding energies deter- 
mine moreover the low transition temperatures. In con- 
trast, on the SC side, the EI is a BEC of tightly bound 
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FIG. 4: Phase boundaries for an EI as a function of a = 
mi/mj. For a<l, the EI is strongly suppressed on the SM 
side because the finite temperature chemical potentials for 
CB electrons and VB holes are different for asymmetric band 
masses leading to breaking of Cooper-type electron-hole pairs. 
On the SC side, on the other hand, transition temperatures 
are only moderately reduced because of the 1/M dependence 
typical for BEC. 



excitons. The higher transition temperatures T c (E g ) on 
this side are however not determined by the larger bind- 
ing energy but by the temperature for which the q = 
exciton state becomes macroscopically occupied (Thou- 
less criterion <|36|) : solid circles). The binding energies 
per se set only the scale for the SC-SM transition 21]. 

The SC-SM transition is driven by screening and Pauli 
blocking due to free charge carriers which, because of 
thermal excitation across a small energy gap, are not only 
available on the SM but also on the SC side. Neglecting 
on the SC side thermally excited charge carriers would 
lead to a much steeper phase boundary (dotted line in 
Fig. [3J and to a T c (0) « 0.45, accidentally identical to 
the guess given for Sr in Ref. || • The relative importance 
of phase space filling vs. screening can be estimated by 
comparing in Fig. [3] the SC-SM boundary (thick dashed 
line) with the hypothetical boundary (thin dashed line) 
arising when only screening is taken into account. In that 
case, the SC-SM boundary is given by B = 0, which is 
equivalent to the ordinary Mott criterion ax • Q s = 1 (see 
Eq. (|32|l ). Contrasting the thin and thick dashed lines 
demonstrates that Pauli blocking cannot be ignored. 

Phase boundaries for asymmetric band masses (a ^ 1) 
are shown in Fig. The differences of the condensation 
mechanisms on the SC and the SM side, respectively, can 
be most clearly seen in this figure, because the transition 
temperature of a BCS condensate of electron-hole pairs 
strongly depends on a whereas a BEC condensate of exci- 
tons does not. In contrast to equal band masses (a = 1), 
where the chemical potentials for the CB electrons and 
VB holes are pinned to —E g /2, as can be seen by inspec- 
tion from Eqs. I|20|l and (|22|l . asymmetric band masses 
lead to chemical potentials which are different for finite 
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FIG. 5: Bound state fraction of the exciton matter above an 
EI with q = 0.0125 as a function of the energy gap and for 
various temperatures. The interceptions with the dotted line 
indicate the energy gaps for which the bound state fraction is 
for the given temperature 50%. 



FIG. 6: Phase boundary of an EI with a = 0.0125 and con- 
tours where the bound state fraction of its halo is 99.9%, 96%, 
70%, and 50%, respectively. The temperature-dependent po- 
sition of the maximum of the bound state fraction is shown 
by the thin solid line. Thin dashed lines indicate the temper- 
atures for which in Fig.|7|the densities are plotted. 



temperatures. Band asymmetry has thus a pair breaking 
effect on Cooper- type electron- hole pairs, similar to the 
effect a magnetic field has on Cooper pairs in supercon- 
ductors. Accordingly, it leads to a strong suppression of 
T c (E g ) on the SM side. On the SC side, on the other 
hand, the EI is supported by tightly bound excitons and 
the 1 /M dependence of the BEC transition temperature 
leads only to a moderate a dependence of T c (E g ). Thus, 
at a ^ 1 the calculated phase boundary deviates sub- 
stantially from the qualitative guess given in Fig. ^ 

Our results seem to suggest that a <C 1 completely 
destroys the EI on the SM side even at T = 0. This is 
however an artifact of our numerics which reaches here 
its limits of accuracy. The EI should be stable at T = 
(but not T 0) for arbitrary mass ratios a. Only 
anisotropies |43| and/or ( mult i-) valleys in the band dis- 
persions separated by 2w ^ K with K a reciprocal lattice 
vector lead at T = to a suppression of the EI on the 
SM side for energy overlaps larger then a critical value. 
Note, multi-valleys in a band which are connected by a 
reciprocal lattice vector K, i.e. identical multi-valleys, 
induce an "effective" mass asymmetry and thus the 
same pair breaking effect for finite temperatures as the 
real mass asymmetry considered above. 

For a <C 1 the EI is almost entirely located on the 
SC side of the SC-SM boundary and within our approx- 
imations a BEC of excitons. Ignoring biexcitons, which, 
for strong mass asymmetry and attractive interactions, 
could interfere on the SC side with the EI 01 j two tem- 
peratures are required to characterize the EI and its exci- 
ton environment (halo of an EI): Tm{E 9 ) where exciton 
formation sets in and T c (E g ) where condensation (phase 
coherence) is finally reached. Above T c (E g ) but below 
Tfy[(Eg), a mixture of excitons, free electrons, and free 
holes exists, the composition of which depends on T and 



Eg, as can be seen in Fig. EI where we plot bound state 
fractions 



7 



(39) 



above T c (E g ) for an EI with a = 0.0125. Here, n{ de- 
notes the part of the total density corresponding to free 
electrons and n\ the part bound in electron-hole pairs 
(excitons) as given by Eqs. (|A.50J) and l|A.51jl . respec- 
tively. 

The bound state fraction indicates on which side of the 
chemical equilibrium e + h ^ X the system is accord- 
ing to the mass-action law 41]. Deep on the SC side, 
the chemical equilibrium is practically for all tempera- 
tures on the side of unbound electrons and holes and the 
bound state fraction is zero. Decreasing the energy gap 
leads to a shift of the chemical equilibrium to the exciton 
side signalled by an increasing bound state fraction which 
assumes a temperature dependent maximum (which can 
be very close to unity) before it abruptly decreases again 
to zero at and beyond the SC-SM transition. The depth 
and width of the maximum strongly depends on temper- 
ature. At low temperatures and large energy gaps the 
bound state fraction acquires a step-like shape (already 
visible by the solid line in Fig. EJ. In that parameter 
range, the exponential tails of the Fermi and Bose func- 
tions contribute to Eqs. I|A.50|) and l|A.51Jl . Moreover, 
B — > B — » R x , Eg — > Eg, and fix w ^ 2 ~ —E g /2, lead- 
ing to 
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1 



1 + exp \J3(R X - Eg/2)} 



(40) 



Thus, as a consequence of the mass-action-law, the bound 
state fraction 7 has a step at E g = 2 ■ Rx for j3E g , /3(E g — 
Rx) 1. 
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FIG. 7: Bound (solid line) and free parts (dashed line) of the 
CB electron density as a function of the energy gap and for 
two temperatures above the EI (dotted lines in Fig.|SJ. Notice 
the very sharp decrease of the bound part of the density when 
the SC-SM transition is approached. 



The phase boundary for the EI with a = 0.0125 to- 
gether with contours of the bound state fraction of its 
halo are displayed in Fig. Almost the whole EI sits 
on the SC side, where it is surrounded by a large exciton 
rich region with bound state fractions significantly above 
50%. Remarkably, the bound state fraction approaches 
almost 100% before phase coherence is established. The 
exciton density in that region is not yet large enough for 
condensation to occur. Notice also, as a consequence of 
Eq. (|40|l . all contours in Fig. [5] start at E g = 2 (exciton 
units). 

Figure shows the exciton dominated part of the 
halo (7 > 0.5). In that region, we expect electri- 
cal resistivity anomalies similar to the ones observed 
in TmSeo.45Teo.55 HI [H. Above T c (E g ) but below 
Tm(E 9 ), for instance at T = 0.071, diminishing E g (in- 
creasing pressure) pushes the system from the free state 
dominated to the bound state dominated regime. In the 
former, the resistivity is expected to decrease with de- 
creasing Eg, because decreasing E g leads to an increase 
of free electrons and holes. The resistivity is here deter- 
mined by the scattering of free charge carriers on imper- 
fections and phonons. In the latter, however, decreasing 
Eg leads not only to an increase of the free charge carri- 
ers but also to a rather strong increase of bound states, 
as can be seen in Fig. Depending on temperature and 
Eg, the bound part of the density can be orders of mag- 
nitude higher then the density of the free charge carriers 
responsible for charge transport. As a result, an addi- 
tional scattering channel is now available, free-bound 
state scattering, and strongly increases the resistivity. 
After a temperature dependent critical E g , the bound 
part of the density decreases with decreasing energy gap. 
Hence, free-bound state scattering diminishes and the 
resistivity decreases again until it changes abruptly to a 



lower value at the SC-SM transition. 

The results we presented in this section are for a 
generic Wannier-type EI. Contact with particular ma- 
terials can be made by specifying the mass ratio a, the 
exciton Rydberg Rx, and the exciton Bohr radius ax 
(or, alternatively, one of the two effective band masses 
mi). 

IV. COMPARISON WITH EXPERIMENT 

We now discuss from the vantage point of our theory 
the experimental data for the pressure-induced SC-SM 
transition in TmSep 45 Te n .55 [12l llrJ. ll4L . Because of 
the fee crystal structure |l2j, the VB has a single max- 
imum at the r point while the CB has minima at the 
three X points of the Brillouin zone. We can formally 
account for the three identical valleys of the CB as well 
as for the spin degeneracy by introducing on the rhs of 
Eqs. fl%|l and (fTT)|l multiplicity factors g\ = 6 and 32 = 2, 
respectively. The remaining equations are unchanged. 

The mass ratio estimates for TmSeo.45Teo.55 range 
from a = 0.01 to a = 0.02 ^5| and optical measurements 
reveal an energy gap of 135 meV and an exciton binding 
energy of 75 meV [l2| . The binding energy is thus rather 
large, which, in addition to the high exciton densities 
(see below), indicates that the effective- mass, Wannier- 
type model may not be quite adequate. Hybridization 
effects characteristic for a mixed-valence material such 
as TmSeo.45Teo.55 may also not be completely captured 
by a two-band model. For our model to mimic the rel- 
evant parts of the band structure of TmSeo.45Teo.55 we 
have to assume that during the excitonic instability, hy- 
bridization does not change much and the composition 
of the two bands is more or less fixed. 

Although the multi-valley CB should favor a metal- 
lic electron- hole liquid the experimentally observed 
increase of the resistivity cannot be explained by an 
emerging metallic phase. We consider therefore the static 
screening model, which a priori favors excitonic phases, 
as an effective model for TmSeo.45Teo.55. The investiga- 
tion of the competition between the two phases should 
be based on an improved description of screening (see 
the discussion following Eqs. (|10|l - H12(l ). taking however 
additional scattering processes into account, which most 
probably stabilize the EI in TmSeo.45Teo.55 against an 
electron-hole liquid. Exciton-phonon scattering, for in- 
stance, could be such a process because the narrow VB 
in TmSeo.45Teo.55 makes the exciton very susceptible to 
phonon dressing. Indeed, phonon signatures have been 
experimentally seen ^5} , but it is not clear whether they 
are the driving force for the excitonic instability or only 
triggered by it. In this respect it should be also men- 
tioned that a static distortion associated with a density 
wave (which would be expected from the indirect gap) 
has not been found experimentally. 



11 




150 _ 



300 



Energy Gap [meV] 

FIG. 8: C omp arison of the experimental data for 
TmSeo.45Teo.55 [UlUI with the phase boundary and the 50% 
contour calculated for an effective model applicable to this 
substance (see text for a description). Open circles, squares, 
and triangles denote the onset of the resistivity anomaly, the 
maximum of the resistivity, and the SC-SM transition, re- 
spectively [T^ . The full diamond on the dotted line indi- 
cates the point where the linear increase of the diffusivity has 
been found in Ref. [l5|. The inset shows the bound part of 
the density, n\, along the thick solid line, which denotes the 
(E g ,T) points for which nj is maximal. At the full circle 
n\ w 1.3 x 10 21 cm" 3 . 

With these restrictions in mind, we present in Fig. [3] 
the phase boundary for a — 0.015, calculated with the 
multiplicity factors given above and scaled with the exci- 
ton Rydberg applicable to TmSeo.45Teo.55, and compare 
it with the phase boundary for the "excitonic insulator" 
given by Wachter and coworkers in Ref. Recall the 
phase boundary was constructed from the anomalous be- 
havior of the electrical resistivity which, in a particular 
pressure and temperature range, first increases, peaks at 
a certain pressure, and then decreases with pressure until 
it terminates in a discontinuous jump. The on-set of the 
increase and the jump have been used to determine the 
entrance points for the "excitonic insulator" from the SC 
(open circles) and the SM from the "excitonic insulator" 
(open triangles), respectively. The E g — T values where 
the resistivity is maximal are given by the open squares. 
Fig. [H] also shows the data point (full diamond on the 
dotted line) where the linear increase of the thermal dif- 
fusivity has been observed. Wachter and coworkers claim 
that at this point TmSeo.45Teo.55 is in a superfluid exci- 
ton state [la ]. 

Apparently, Wachter and coworkers ^5| distinguish be- 
tween an "excitonic insulator" at high temperatures and 
a superfluid exciton state at low temperatures which, 
however, is misleading, because the term "excitonic insu- 
lator" has been originally introduced to describe, on the 
SM side, a superfluid electron-hole condensate, similar 
to a superfluid condensate of Cooper pairs in a super- 
conductor, and, on the SC side, an exciton condensate, 
similar to superfluid He 4 H 0, 0, S E3 • Although the 



precise meaning of the terms "excitonic insulator" and 
"exciton condensation" changed in the course or their 
work and, in fact, was never fully in accordance with the 
original definitions, in more recent publications [LH [l5| 
they denote by the former an exciton liquid and by the 
latter a gas/liquid transition. The transition to the su- 
perfluid state they consider to be similar to the one in 
liquid He 4 [la]. 

Based on the results presented in the previous section, 
we suspect that the resistivity anomaly in TmSeo.45Teo.55 
is due to free-bound state scattering in the exciton rich 
part of the halo of an EI whereas the linear increase of 
the thermal diffusivity could signal an EI, understood 
as a macroscopic, phase coherent condensate of excitons. 
If our hypotheses are correct, we should be able to (i) 
relate the phase boundary of the "excitonic insulator" to 
the above 50% contours of the bound state fractions and 
(ii) find the experimental signature for superfluidity at a 
temperature where we would expect the EI. 

The experimentally determined exciton Rydberg for 
TmSeo.45Teo.55 puts the exciton rich part of the halo in 
the temperature range T ~ 50 - 180 K, close to the tem- 
perature range where the resistivity anomaly and thus 
the "excitonic insulator" has been observed (open trian- 
gles, squares, and circles in Fig. [SJ. The bound part of 
the density in this range has also the correct order of 
magnitude. At T — 53 K and E g = 64 meV (full circle in 
Fig. |HJ), for instance, we find, using mi = 2.1m e for the 
CB mass (m e bare electron mass), n\ « 1.3 x 10 21 cm' 3 
which compares favorably with the experimentally esti- 
mated exciton density of 3.9 x 10 21 cm -3 0]. The bound 
part of the density reaches its maximum along the thick 
solid line in Fig. [S] As can be seen in the inset, below 
100 K, i.e. within the bound state dominated region, 
n\ increases sharply. An improved theory, taking inter- 
actions between electron- hole pairs into account, could, 
when the repulsive part of the interaction dominates, dis- 
play a gas/liquid transition in that temperature range. 
Hence, an exciton fluid, as proposed in Ref. 0], is con- 
ceivable below 100 K. 

Below 50 K we find the EI, that is, the macroscopic, 
phase coherent condensate, very close to the 20 K (dot- 
ted line) where the linear increase of the diffusivity was 
observed at 13 kbar. From the associated isobar, we esti- 
mate the energy gap for that data point to be « 20 meV 
(full diamond). Wachter and coworkers' claim that this 
increase is due to the second sound of a superfluid is sup- 
ported by our theory because the EI is entirely on the 
SC side. Hence, it constitutes a BEC of bosonic excitons, 
and not a BCS condensate of loosely bound electron-hole 
pairs, which may or may not display superfluidity pjif . 
Without detailed knowledge of the interaction between 
electron-hole pairs, we cannot calculate superfluid prop- 
erties. Nevertheless, the fact that the EI is on the BEC 
side makes it plausible that an improved theory would 
yield superfluidity, provided the interaction between ex- 
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citons is dominantly repulsive. 

The temperature dependence of the bound state frac- 
tions is consistent with values obtained form ab-initio 
numerical simulations of mass asymmetric electron-hole 
plasmas |45j . The approximate T-matrix H26JI on which 
the calculation of the bound state fractions is based cap- 
tures therefore the essential physics. The interpretation 
of the experimental data is however hampered by two 
main difficulties. 

First, the effective-mass model and the resulting 
Wannier-type EI are, strictly speaking, not applicable to 
TmSeo.45Teo.55. The exciton binding energy is too large 
compared to the energy gap. As a result, the bound state 
dominated region in Fig.|Hlcomplctely exhausts the maxi- 
mal energy gap of 135 meV available at ambient pressure. 
Hence, the free state dominated region, where we expect 
a decrease of the resistivity with decreasing energy gap, 
occurs at energy gaps which cannot be realized. 

Second, the bare energy gaps E g , which we employ 
as a control parameter, are not necessarily the energy 
gaps given by Wachter and coworkers. Using the linear 
relation E g (p, T = 300 K) = E g (p = p a ,T = 300 K) + 
(dE g /dp)(T = 300 K)(p — p a ), with p a the ambient pres- 
sure, they convert (pressure, temperature) points, which 
they control experimentally, to (energy gap, tempera- 
ture) points assuming the energy gap to be tempera- 
ture independent and the closing rate to be temperature 
and pressure independent |l2l fl3j |. At finite tempera- 
tures, however, bands are partially filled. The band gap 
renormalization to which the partial filling gives rise to 
leads not only to a strong temperature dependence of 
the energy gap but most probably also to a complicated 
pressure and temperature dependence of the closing rate. 
These effects were not incorporated in the assignment 
(pressure, temperature) <-> (energy gap, temperature) 
and perhaps the reason for the slight 'left-turn' of the 
experimental curves. 

V. CONCLUSIONS 

We adopted a BEC-BCS crossover scenario to cal- 
culate the phase boundary of a pressure-induced EI 
using an isotropic, effective-mass model to describe 
the two bands involved in the excitonic instability 
and the dominant-term-approximation for the (statically 
screened) Coulomb interaction. 

Assuming pressure to affect only the energy gap, we 
derived, within a meanfield approximation, a linearized 
gap equation for the anomalous selfenergy (order param- 
eter) and determined, for a given energy gap E g , the 
transition temperature T c (E g ) as the highest tempera- 
ture for which a non-trivial solution exists. The kernel 
of the gap equation depends on the chemical potentials 
for CB electrons and VB holes and on the screening wave 
number, all three have been determined self-consistently 



as a function of T and E g . 

We pointed out that because pressure directly controls 
the chemical potentials, which moreover are constrained 
by charge neutrality, the meanfield theory is even for fi- 
nite temperatures on par with T-matrix based crossover 
theories which include pair states in the equation for the 
chemical potential. We demonstrated by physical reason- 
ing and calculation that in our case, pair states do not 
affect the chemical potentials as long as one stays within 
the framework of noninteracting pairs. As a result, the 
transition temperatures obtained from our meanfield the- 
ory coincide on the SC side with the transition tempera- 
tures for a BEC of non-interacting bosonic excitons. As a 
consistency check, we also determined T c {E g ) on the SC 
side from the normal phase electron-hole T-matrix using 
the Thouless criterion. 

Our main results are as follows: For equal band masses 
we found a very asymmetric phase boundary T c (E g ), 
with a steeple on the SC side and a long tail on the 
SM side, directly reflecting the different character of the 
EI on the SC and SM side, respectively. Whereas on the 
SM side the EI comprises a BCS condensate of loosely 
bound Cooper pairs, the EI is a BEC of tightly bound 
excitons on the SC side. On the SM side, no simple ex- 
ponential dependence of the transition temperature on 
the band overlap \E g \ exists, as perhaps expected from 
the formal analogy to the BCS theory of superconduc- 
tivity. We found mass asymmetry to strongly suppress 
pairing on the SM side because the different tempera- 
ture dependencies of the chemical potentials for the CB 
electrons and VB holes have a pair-breaking effect on 
Cooper-type electron-hole pairs. The transition temper- 
atures on the SC side, on the other hand, are only mod- 
erately modified, because the center-of-mass motion of 
excitons, which determines the transition temperatures 
on that side, depends only algebraically on the total mass 
(T c ~ M ~ 1 ). For very asymmetric band masses the EI is 
therefore rather fragile on the SM side and for finite tem- 
peratures almost entirely located on the SC side of the 
SC-SM transition. On the SC side, we also investigated 
the exciton density (more precisely, the part of the total 
density bound in excitons) in the region above T c which 
we called halo. Because of the mass-action-law, the exci- 
ton density strongly depends on temperature and energy 
gap. As a consequence, we expect pronounced resistivity 
anomalies in the halo due to free-bound state scattering. 

Within the limitations of our model, we attempted 
an analysis of the experimental data for the pressure- 
induced SC-SM transition in TmSeo.45Teo.55 El 
Hil llij . From our theoretical point of view, the phase 
boundary, which was constructed from electrical resistiv- 
ity data, traces most probably the exciton rich surround- 
ings of an EI, where excitons dominate the total density, 
provide abundant scattering centers for the charge carri- 
ers, and induce the observed resistivity anomaly. Since 
we neglect interactions between excitons, we cannot de- 
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cide whether a gas/liquid transition is also involved in 
the resistivity anomaly ^5|. The linear increase of the 
thermal diffusivity which was taken as a signal for super- 
fluidity occurs at a temperature where we would expect 
the EI. Most importantly, however, the EI and the data 
point are on the SC side of the SC-SM transition, where 
the EI, that is, the macroscopic, phase coherent exciton 
state, consists of tightly bound bosonic excitons, which 
should support superfluidity when the repulsive core of 
the exciton-exciton interaction prevails over the attrac- 
tive tail. T hus, the long journey of Wachter and cowork- 
ers 0, 0, Q, 0| may indeed have culminated in 
the first observation of an exciton condensate. 

We could not unambiguously proof the existence of a 
superfluid exciton condensate in TmSeo.45Teo.55, most 
notably, because we did not allow the EI to compete 
with a metallic electron-hole liquid and because we did 
not account for the possibility of density waves. Al- 
though neither is supported by experimental data, both 
are in principle possible in a multi-valley band struc- 
ture. Thus, to further substantiate the experimentalists 
claim from the theoretical side, a unified treatment of 
excitonic and electron-hole liquid phases, accounting also 
for the possibility of density waves, is clearly necessary, 
ideally for a model which avoids the effective-mass ap- 
proximation, includes exciton-phonon scattering (which 
could stabilize excitonic phases against an electron-hole 
liquid), and an approximation scheme which takes the 
interaction between electron-hole pairs (or excitons) into 
account. The latter is important for the calculation 
of transport coefficients, which could be directly com- 
pared with experimental data, as well as for clarifying 
the role of biexciton formation, which, for large mass 
asymmetry, competes with the formation of the EI. Be- 
sides diagrammatic techniques |2n ] and systematic mode- 
couplin g appro aches (23i l29| , numerical simulation tech- 
niques |3JJ, |45| could proof very useful in this respect. 
In any case, the comparison between theory and experi- 
ment depends on the assignment (pressure, temperature) 
<-> (energy gap, temperature), which is rather subtle. 
To study the SC-SM transition along 'iso-Eg' instead 
of isobars would eliminate the uncertainties associated 
with this critical assignment and thus open the door 
for a more rigorous analysis of the excitonic phases in 



TmSe .45Te ., 



APPENDIX A: CB ELECTRON AND VB HOLE 
DENSITIES IN T-MATRIX APPROXIMATION 



In this appendix we determine the CB electron and 
VB hole densities as a function of the bare energy gap 
and the temperature. Thereby, we verify that even in the 
T-matrix approximation the chemical potentials for CB 
electrons and VB holes are, in leading order, not affected 
by exciton states. 

Our starting point are the well-known formulae 



n-2 



dk du 

(2?r) 3 2tt 

dk du 

(2tt) 3 2^ 



A u (k,w)n P H, (A.41) 
A 22 (k,-uj)n F (uj) (A.42) 



for the CB electron and VB hole density, respectively. 
The spectral functions An are calculated with the self- 
energies shown in Figs. an d (El 3 )- Expanding the 
spectral functions with respect to Tu = — 2ImSa leads 

to mi 

Aii(k,w) = 27rZ«(k)$(w - e^k)) 

+ r^k,^)— ^— T (A.43) 



ctej(k) lu - ej(k) 



with the quasi-particle weight 
Z«(k) = 1 



de^k) J 2™-e 4 (k) ^ ' 

and the renormalized single particle dispersion 

ei(k) = ei (k) + (k) + ReE£(k, e,(k) + i v ). (A.45) 

Let us first consider the CB electron density. Assuming 
the internal propagators in the selfenergy to describe free 
quasi-particles with a renormalized dispersion gj but a 
quasi-particle weight unity, the selfenergy arising from 
the T-matrix reads 



Sfi(k,w„) = / ~ V) A12 (k, k; q, iQ. n )Q 22 (k - q, iu n - iSl n ) 
J ^ 13 ml 



^ fr(|k - £q|)] ^ n l iE ^ q)) + t ? {e t~ ^ " ■ ( A ^) 



(2tt) 3L ^' Ex(q)-s 2 (k-q)+tJ, 2 -iu n 

where we introduced the Bose function ub and the exciton dispersion Ex(<l) — E g — B + q 2 /2M. To obtain Eq. (|A.46|I . 
we approximated the renormalized dispersions £j, that is, the solutions of Eq. I|A.45|I . by Eq. I|14|) . For the electron-hole 
T-matrix we used Eq. (12(it with Eq. I|34|l for the renormalized exciton propagator. Thus, after analytical continuation 
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to real frequencies, we obtain 



rn(k, w) = J [,g(|k - ^q|)] 2 ^x [n B (E x (q)) + n F (e 2 (k - q) - ^ 2 )]S(E x (q) - e 2 (k - q) + ,i 2 - w)(A.47) 



for the imaginary part of the selfenergy on the real fre- 
quency axis and, using a Kramers-Kronig relation, 



with 



ReEfi(k,u;) = P J 



duj' r n (k,u}') 

2lt LU — LU' 



(A.48) 



for the real part of the selfenergy on the real frequency 
axis. 

If we now insert Eqs. (|A.47j) and (|A.48jl into Eq. (|A.43ll . 

recall once again Eq. (|14[) . and perform the k-integration 
in Eq. (|A.41() . we find for the CB electron density 



with a term arising from free electrons (in the sense de- 
fined above) 



dk 



n F (ei(k) - m) 



(A.50) 



and a term 



dpdq 



D 



p 

2m 



B 



P2 
2m 



X 2 (p)Z x F 1 (p,q) (A.51) 



due to the coupling of CB electrons and VB holes. The 
statistical factor is given by 

Fi(p,q) = n B (-Sx(q))[l-n F (ei(p+ -rjq) -Mi) 



n F (£2(p 



m 2 



ft , mi ^ 
n F (ei(p+ -^-q) 



x n F (e 2 (p 



m 2 



M2)] 

Mi) 
Ma)< 



(A.52) 



To derive Eq. I|A.51(1 we transformed the k-integration in 
Eq. I|A.41(I to an integration over the relative momentum 
p = k (mx/M)q when Tu was part of the integrand 
and used definition (|28|1 of the form factor. 
In the limit f3E g <C 1, Eq. (|A.51|) reduces to 



J^MExiq)) 



(A.53) 



because Fi(p,q) — > nB(Ex(q)), Z x — > 1, and the inte- 
gration over p yields unity. The latter, a consequence of 
the optical theorem which the T-matrix has to obey, can 
be only ensured within a skeleton expansion. 

The VB hole density can be calculated in a similar way. 
The result is 



n 2 



-J 



(A.54) 



-J - 



dk 



■7lF(e2(k) - (Mi) 



(A.55) 



and a bound part n\ which is identical to n\ given in 
Eq. (|A.51J| . This is not unexpected, because for each CB 
electron participating in a bound state, a VB hole has to 
do so too. 

As a consequence, the charge neutrality condition, 
m = fi 2 , which fixes the chemical potentials for CB elec- 
trons and VB holes, reduces to n{ = n^. Neglecting 
the modifications of the single particle dispersions due 
to excitons, that is, ignoring interactions between exci- 
tons, which are higher order effects, the charge neutrality 
condition in the T-matrix approximation reduces to the 
mcanfield charge neutrality condition. 

With the chemical potentials determined from 
Eqs. J20I) and Eqs. (|A.50|) and (|A.51|) can be used 
to calculate, above the critical temperature T c (E g ), the 
free and bound parts of the CB electron density and thus 
the bound state fraction ||3"§1) . 



APPENDIX B: SEPARABLE APPROXIMATION 
FOR THE NORMAL PHASE ELECTRON-HOLE 
T-MATRIX 



Here we derive a separable approximation for the nor- 
mal phase electron-hole T-matrix Ai2(k, k'; q, ifl n ). In 
the vicinity of the exciton state the T-matrix factorizes 
with respect to k and k'. Thus, the separable approxi- 
mation is a good approximation for energies not too far 
from the exciton energy. Since this energy range deter- 
mines the physics of an EI, the separable approximation 
should be sufficient for our purpose. 

The approximation is based on an expansion of the nor- 
mal phase electron-hole T-matrix in terms of the eigen- 
functions of an auxiliary Lippmann-Schwinger equation 
which takes into account screening due to spectator par- 
ticles not participating in the bound state but ignores 
the spectators' Pauli blocking. The phase space filling 
responsible for Pauli blocking is incorporated in a second 
step through a selfenergy for the exciton propagator. To 
avoid the numerical solution of the auxiliary Lippmann- 
Schwinger equation we approximate, as far as the calcula- 
tion of the T-matrix is concerned, the screened Coulomb 
potential by the Hulthcn potential. Eigcnfunctions and 
binding energies are then analytically available. 
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The starting point is Eq. (|24|) . In a first step, we break up the screened ladder equation into two equations. In an 
obvious notation 

A^[Q-Q^}A 12 , 



A 12 = Ag> 



A 



(0) 



= V S + V S Q^4 



(<>) 

2 



with 



Q<°>(k;q,ifi n ) = 



(A.56) 
(A.57) 

(A.58) 



Jg -r £i(k) +e 2 (k - q) - ' 

Note, this is an exact rewriting of Eq. H24|) . with the advantage that screening and Pauli blocking are now treated 
separately. Equation ljA.571) introduces an auxiliary vertex which considers screening but not Pauli blocking. The full 
vertex is then obtained from Eq. I|A.56[1 . 

Introducing relative momenta p according to k = p + (mi/M)q, Eq. i|A.57|) becomes 

771l i . 

-q;q,iS2 n ) = V s {p-p ) 



Ala (p + ^-q,p 



M 



dp" 



V S (P- 



(2^) 3 E a 



2M 



. (o) . ., m\ , 
A ,' p" H -q, p' 



777l 
~M 



which, with the identification 



A - (P+ M q ' P + M 



reduces to 



q; q, «^«) = T(p, p', m n - E g 



T(p,p';ifi n ) = V.(p-p') + J |Pl^-^T(p".p' : /o„.. 



2 AT 



q;q,ifi n ), (A.59) 



(A.60) 



(A.61) 



2m 



Equation l|A.61|) is the T-matrix equation in the center-of-mass frame of an electron-hole pair interacting via a screened 
Coulomb potential. The medium occurs here only through the screening wave number q s . 

To solve Eq. (|A.61() we employ the unitary pole expansion for T(p, p'; ifl n ), using eigenfunctions of the corresponding 
Lippmann-Schwinger equation at fixed energy as a basis 46]. In particular, we take the eigenfunctions at energy 
iQ n — > —B, with —B the energy of the lowest bound state. Truncating the expansion after the first term yields (for 
an isotropic system) 

T(p,p';in n ) = g(p) ■ d x (iQ n ) ■ g(p') (A.62) 

with 

-1 



d x (z) = 



(A.63) 



z + B 1 

which can be interpreted as a screened exciton propagator, and form factors defined in Eq. I|28|). For the Hulthcn 
potential the screened exciton wavefunction x{p)> an d thus the from factor g(p), as well as the screened binding energy 
B can be obtained analytically and are given by Eqs. Q31f) and Q.'j2H . respectively. 
Combining now Eq. l|A~62|) with Eq. (|A~60|) leads to 

(A.64) 



A^(P + ^q, P' + ^jq q «1„) = 9 {p) • d X (iSl n - E g - ■ g(p') 



for the auxiliary vertex. To obtain the full vertex, we write A12 in a similar form, 
Ai 2 (p+ -Tjq,p' + -tj q;q,«n„) = g(p) ■ D x (q,i£l n - E g 



2M 



(A.65) 



but with a renormalized exciton propagator Dx{q, ifl n ) instead of the screened exciton propagator dx(iQ n )- Inserting 
Eqs. (|A.64I) and (|A.65|) into Eq. (|A.56(I . leads to a Dyson equation, 

£> x (q,if2 n ) = dx(itt n ) +d x (ifl n )M x {ci,in n )Dx(ci,in n ), (A.66) 

for the renormalized exciton propagator Dx(q, iQ n ) with the selfenergy Mj(q, zf2 n ) given in Eq. (|29p. The selfenergy 
takes Pauli blocking into account. Switching in Eq. (|A.65|) from the relative momenta p and p' back to the original 
momenta k and k', we recover Eq. I|26[l for the full vertex. 
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